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Ab initio molecular-dynamics simulations have been used to investigate the structure, dynamics 
and electronic properties of the liquid alloy Agi-^iSe^, at 1350 K and at the three compositions 
X — 0.33, 0.42 and 0.65. To provide a point of reference, calculations are also presented for the 
equilibrium structure and the electronic structure of the a-Ag2Se crystal. The calculations are 
based on density-functional theory in the local density approximation and on the pseudopotential 
plane-wave method. For the solid, we find excellent agreement with experiment for the equilibrium 
lattice parameters and the atomic coordinates of the 12-atom orthorhombic unit cell, and we present 
an analysis of the electronic density of states and density distribution. The reliability of the liquid 
simulations is confirmed by detailed comparisons with very recent neutron diffraction results for the 
partial structure factors and radial distribution functions (RDF) of the stoichiometric liquid Ag2Se. 
Comparison with the predictions of an empirical interaction model due to Rino et al. is also given 
for ^-Ag2Se. The ab initio simulations show a dramatic change of the Se-Se RDF with increasing 
Se content. This change is due to the formation of Se clusters bound by covalent bonds, the Se- 
Se bond length being almost the same as in pure c-Se and ^-Se. The clusters are predominantly 
chain-like, but for higher x there is a significant fraction of 3-fold coordinated Se atoms. It is shown 
that the equilibrium fractions of Se present as isolated atoms and in clusters can be understood 
on a simple charge-balance model based on an ionic interpretation. The Ag diffusion coefficient 
in the simulated stoichiometric liquid is consistent with experimental values measured in the high- 
temperature superionic solid. The Ag and Se diffusion coefficients both increase with Se content, in 
spite of the Se clustering. An analysis of the Se-Se bond dynamics reveals surprisingly short bond 
lifetimes of less than 1 ps. The electronic density of states (DOS) for ^-Ag2Se strongly resembles 
that of the solid. Some of the changes of DOS with composition arise directly from the formation of 
Se-Se covalent bonds. Results for the electronic conductivity a obtained using the Kubo-Greenwood 
approximation are in adequate agreement with experiment for £-Ag2Se, but for the high Se contents 
the simulation results for cr are 3-4 times greater than experimental values. Possible reasons for this 
are discussed. 
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I. INTRODUCTION 

The study of hquid binary aUoys has been an ex- 
tremely fruitful source of insights into the relations be- 
tween structure and electronic properties in condensed 
matter. There has been a vast amount of experimental 
work on metallic and semiconducting binary liquids, and 
it has been well known for many years that tteir proper- 
ties often vary dramatically with compositjonEllj. Famous 
cases are the Cs-Au and Mg-Bi alloysoQ, in which the 
pure elements are good metals, but nevertheless the sto- 
ichiometric mixtures (CsAu and Mg3Bi2) have very low 
conductivities, and their structures are characteristic of 
molten salts. These effects arise from the electronega- 
tivity difference between the elements, and the resulting 
charge transfer, partial ionicity and heterocoordination. 
In systems where one of the elements is a semiconductor 
in the liquid state - for example alloys of metals with S or 
Se - even richer behavior can be expected, since variation 
of composition should change the bonding from metallic 
through partially ionic to covalent. We report here a set 
of simulations of the liquid Ag-Se system performed us- 
ing ab initio molecular dynamics (AIMD), which we have 
used to explore these effects. 

The electrical conductivity a of the Agi_a;Sea; system 
has been measured over most of the jcomposition range 
at temperatures from 973 to 1573 Kdil. As the com- 
position goes from pure £-Ag to the stoichiometric al- 
loy Ag2Se, (T decreases from the typically metallic value 
of ~ 50000 f2~^cm~^ to typical semiconducting values 
of ~ 500 fi^^cm^^. This strongly suggests the forma- 
tion of a pseudogap in the electronic density of states 
at the stoichiometric composition, with the Fermi energy 
lying in this gap. As the Se content is further increased, 
a varies rather little over the range between x = 0.33 
and X = 0.65, and then descends to very low values 
of ^ 1 il~^cm~^ at 700 K as the Ag content goes to 
zeroEI. The propiecties of pure ^-Se have been very ex- 
tensively studiedE3, and it is known to have a gap of 
ca. 1.9 eV in the density of states just above its melting 
point (~490 K), so that it is a wide-gap semiconductor. 
This gap decreases with increasing temperature, and ap- 
pears to be ca. 0.7 eV at 1350 K and 100 barO. 

Pure ^-Ag has the rather closed-packed structure and 
high coordination number (ca. 12) expected from its 
f.c.c. crystal structureEj. Until very recently, the struc- 
ture of £-Ag2Se had been measured ouh at the level 
of the neutron-weighted structure factoJlj, but neutron 
diffraction combined with isotope variation has now been 
used to measure the partial structure factors and hence 
the partial radial distribution functionaij. These re- 
veal the heterocoordination characteristic of ionic or par- 
tially ionic liquids, as will be discussed more fully below. 
Diffraction studies-piiJ-Se have been made by several in- 
dependent groupstjilZl, and it is well established that the 
average coordination of the atoms is close to 2.0 over a 
wide range of temperatures and pressures from the melt- 



ing point up to 1773 K and 815 bar. It is widely believed 
that this coordination indicates the presence of extended 
chain-like struct|i,iy.ps, and there is support for this from 
simulation workOEj. These facts imply that there must 
be major changes of atomic ordering as the composition 
is varied from pure ^-Ag to pure ^-Se. However, nothing 
is yet known about the structure of the alloys except at 
the stoichiometric composition. 

The aim of the work reported here is iiji|_Lise ab initio 
molecular dynamics simulation (AIMD)E3£j to investi- 
gate the structure, the atomic dynamics and the elec- 
tronic structure of the liquid Ag-Se alloys at different 
compositions. The issues we want to address are: the re- 
lationship between the electronic structure of the liquid 
and the solid; the characteristics of the electronic density 
of states and electron distribution as a function of compo- 
sition; the variation of liquid structure with composition, 
and particularly how the apparently ionic structure of £- 
Ag-Se goes over to the chain-like structure of ^-Se; and 
the way this liquid structure is related to the dynamics 
of the ions. At present, there is almost no experimental 
information about any of these questions. Of course, be- 
fore we address any of these issues, we have to show that 
our AIMD techniques are able to give a faithful represen- 
tation of the real material. Because of this, we will begin 
by presenting calculations on the a-Ag2Se crystal, which 
has a rather complex structure and gives a good test of 
our methods. Our study of the crystal also provides an 
essential reference point for discussing the properties of 
the liquid. 

AIMD is ideally suited to this type of problem. The 
energetics of the system and the forces on the ions are 
calculated by first-principles quantum methods, with no 
adjustable parameters, and without the need for the em- 
pirical interionic potentials used in earlier-simulations on 
many liquid mixtures, including £-Ag2SeE3. In addition, 
the thermal-equilibrium structure, atomic dynamics and 
electronic structure are all calculated together in a com- 
pletely unified and self-consistent way. AIMD techniques 
have been used to investigate a number of other liquid 
metals and sexaiconductors, includinjg-the binary alloys 
K-SlB, Na-Snld, Cs-PbB and Ga-SeB. Preliminary re- 
sults of the work-presented here have already been pub- 
lished elsewhereE3. The rest of the paper is organised as 
follows. In Sec. II, we summarize the simulation tech- 
niques, and we then present in Sec. Ill our calculations 
on the equilibrium structure and the electronic structure 
of Q!-Ag2Se. Our simulation results on the liquid Ag-Se 
alloys are reported in Sec. IV, where we examine succes- 
sively the structure, the atomic dynamics and the elec- 
tronic properties, including the density of states and the 
electronic conductivity. We present a brief, discussion of 
the relation between atomic ordering and electronic prop- 
erties in Sec. V, and our conclusions are summarized in 
Sec. VI. 
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II. COMPUTATIONAL METHODS 

Our ab initio molecular dynamics (AIMD) teclipijque is 
closely related to the Car-Parrinello approachoO. We 
use density functionaJ-|tbaory within the local density 
approximation (LDA)c3~e3. Only the valence electrons 
are treated explicitly, the interaction of the valence elec- 
trons with the atomic cores being represented by norm- 
conserving non-local pseudopotentialsEllES. The simula- 
tions are performed in a periodically repeated cell with 
the wave functions expanded in plane waves. The basis 
set includes all plane waves whose kinetic energy is less 
than a chosen cut-off energy i?cut • 

An important difference between our calculations and 
the original Car-Parrinello approach is that instead of 
treating the electronic degrees of freedom by 'fake dy- 
namics' we relax the electrons to the Born-Oppenheimer 
surface for every ionic— coHfiguration, using conjugate- 
gradients minimizatior]E3~Ea. As in the original scheme, 
we calculate the forces acting on the ions via the 
Hellmann-Feynman theorem, and these are then used to 
integrate the classical equation of motion of the ions. It 
is computationally more expensive to bring the electrons 
to the Born-Oppenheimer surface than to do one step 
in the Car-Parrinello method, but the extra cost is com- 
pensated by the fact that we can use a larger molecular 
dynamics time step. To handle the semi-metallic nature 
of the system, we use Fermi-surface smearing, with the 
electronic occupai ti ft n .. n umbers treated as auxiliary dy- 
namical variablesE^lEaO. A detailed report on how this 
is dyy in practice has been presented in our work on 

Although AIMD has become a standard technique, 
very few simulations have been reported on_sizstems con- 
taining transition or post-transition metalsEjO. Because 
of the important role of the d-electrons, it is essential 
to include them explicitly, and this brings two kinds of 
problems. Firstly, the number of occupied orbitals is sub- 
stantially increased. Secondly, the compact nature of the 
d-orbitals means that a large cutoff Ecut is needed, and 
this implies a large plane- wave basis set. Our AIMD 
simulations on the Ag-Se system were performed with a 
version of the CETEP codecd running on a Cray T3D 
parallel supercomputer. We note that the version of the 
code used here reaches the ground state by minimizing 
with respect to all bands simultaneously, rather than by 
the band-by-band method of earlier versions. In the liq- 
uid simulations, we needed about 12 minimization steps 
for every ionic configuration to converge the total free 
energy within 10~^ eV/atom. 

Technical details of our ab initio pseudopotentials are 
as follows. The Ag pseudopoteKtial has been optimised 
using the method of Lin et alsB, iwhich is a refinement 
of the scheme due to Rappe et alxn. The Se pseudopo- 
tential does no|L-|require optimisation and the standard 
Kerker methodo suffices. For Ag, all states below Ad, 
5s and 5p are treated as core states. The s and p com- 



ponents of the Ag pseudopotential were generated using 
the atomic configuration 4:d^'~'5s'^^^^5p'^■^^ , and the d com- 
ponent using the configuration 4d^''5s°'^. The core radii 
were 2.0, 2.0 and 2.5 a.u. for the s, p and d components 
respectively. For Se, we used As^Ap'^ for the s and p waves 
and As^Ap^ '^^Ad^-'^^ for the d wave; the core radii were 
chosen to be 2.0, 2.0, and 2.3 a.u. for the s, p and d com- 
ponents respectively. We work with tlie. pseudopotentials 
in Kleinman-Bylander separable formEJ, with the s-wave 
treated as local; the non-local parts of the pseudopoten- 
tials are treated in real spacec3. We use a plane-wave 
cut-off of 400 eV. As a preliminary test of our pseudopo- 
tentials, we have performed calculations on the low tem- 
perature phase a-Ag2Se, which we present in more detail 
in the next section. We also tested the pseudopotentials 
in our work on AgClEll and the liquid Ga-Se alloyn3. 



III. SOLID AgaSe 

At atmospheric pressure silver selenide, Ag2Se, exists 
in two polymorphic forms. The high temperature form, 
/3-Ag2Se, is stable above 406 K. The structure is b.c.c. 
with spacing a = 4.983 A, and the unit cell contains two 
Ag2Se units. The Se atoms form a stable b.c.c lattice, 
while the Ag atoms are-statistically distributed among 
several interstitial sitealS'CJ. In this structure, all the 
Ag atoms are mobile, and in fact Q;-Ag2Se is a tiqpical 
superionic conductor similar to a-AgI and Q;-Ag2SE£rE3. 
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FIG. 1. Arrangement of Ag and Se atoms in solid a-Ag2Se. 
Ag and Se atoms are represented by small dark spheres and 
large light spheres, respectively. 
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Ag2Se transforms from the /3 phase to a low tempera- 
ture form, a-Ag2Se, at 406 K. The structure of a-Ag2Se 
has been resolved by X-ray powder diffraction as or- 
thorhombic with cell constants a — 4.333, h — 7.062 and 
c = 7.764 A. The space group is_P2i2i2i and there are 
four Ag2Se units in the unit celS. A three dimensional 
representation of the structure is shown is in Fig. |l|. The 
four Se atoms are equivalent by symmetry, and they al- 
most lie on two planes perpendicular to the b axis. There 
are two inequivalent Ag sites. The first type, Agi, lie just 
above or below the planes of the chalcogen atoms. These 
atoms are tetrahedrally coordinated by Se atoms, three 
in the nearest Se plane and one in the next, at distance 
2.62, 2.79, 2.86 and 2.71 A, respectively. The other sil- 
ver atoms, Agii, lie almost halfway between the planes 
of chalcogen atoms. The arrangement of the Se around 
the Agn atoms is almost triangular, at distances 2.72, 
2.74 and 2.81 A, with two second neighbors at 3.28 and 
3.50 A. Each Se atom is surrounded by seven Ag atoms, 
or nine if the second neighbors are taken into account. 

Silver is well known to exist in various solid-state com- 
pounds with coordination numbers ranging from 2 up to 
8Ej. Examples ofj-cpordination 2, 3, and 4 are given by 
Ag20El, AgaAsSaEa and Agll^. In this respect the binary 
compounds of silver with the chalcogen elements, Ag2X 
(X=0, S, Se and Te), are remarkable since the coordina- 
tion of silver varies from linear to-trigonal bipyramidal 
when going down the Via columnE^. With a valence of 
one for silver these compounds may be regarded as ionic 
which makes these low coordinations all the more strik- 
ing, since ionic bonds usually imply higher coordination. 
The stability of the Ag d^^ ion in such compounds may 
be attributed to the mixing of the silver d and s states 
which occurs wh en the atom is placed in a low-symmetry 
environmentLOH. 

The electronic properties of a-Ag2Se are those of a 
semiconductor with a small energy gap. The exact value 
of the energy gap, is however controversial, and_ijalues 
ranging from 0.02 to 0.22 eV have been publishedQO. 

TABLE I. Comparison of calculated and experimental 
crystallographic parameters for a-Ag2Se (experimental val- 
ues in parentheses). The lattice constants are given in A . 
The internal parameters x, y, z are the coordinates of the 
atoms in units of a, b and c. The three other positions equiv- 



alent by symmetry are 


given by ( i 
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Lattice constants 


a 


b 


c 




4.218 


6.949 


7.649 




(4.333) 


(7.062) 


(7.764) 


Internal parameters 


X 


y 


2 


Agi 


0.109 


0.367 


0.455 




(0.107) 


(0.366) 


(0.456) 


Agn 


0.729 


0.027 


0.362 




(0.728) 


(0.029) 


(0.361) 


Se 


0.362 


0.245 


0.154 




(0.358) 


(0.235) 


(0.149) 



The study of a-Ag2Se, with three lattice parameters 
and nine internal parameters to relax, is a stringent test 
of the quality of our pseudopotentials. All the calcu- 
lations on a-Ag2Se were done with a plane-wave cut- 
off of 400 eV. This cut-off yields well-converged proper- 
ties of the fully relaxed structure, since increasing £^cut 
to 600 eV changed the forces on the atoms and the 
stresses on the unit cell only by a negligible amount. 
The Brillouin zone (BZ) sampling was p^formed using 
a Monkhorst-Pack special fc-points meshEa. We also had, 
to introduce a Gaussian smearing of the Fermi surfaceES 
due to the absence of a real gap in the system. We had 
to use a (4x4x4) grid (equivalent to 8 /c-points in the 
irreducible part of the BZ) for the energy per atom to be 
converged within a few meV. The energy broadening we 
used was 0.2 eV. 

We performed a full structural relaxation of the crys- 
tal structure by varying the three lattice parameters and 
the nine independent internal parameters. In Tab. |, 
the calculated lattice constants and internal parameters 
are compared with experiment. The lattice parameters 
are in satisfactory agreement with experiment, being too 
low by ~2%, which is typical of LDA calculations. The 
experimental internal parameters for the Ag atoms are 
very well reproduced. For the Se atoms however, the 
relaxations were larger, and this leads to larger discrep- 
ancies with experiment. The resulting Ag-Se distances 
compared with experimental values are (A units): 2.68 
(2.62), 2.72 (2.78), 2.89 (2.86), 2.80 (2.71) for Agi and 
2.67 (2.74), 2.74 (2.72) and 2.80 (2.81) for Agn. 
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FIG. 2. (a) Electronic density of states of solid Q:-Ag2Se 
and (b) local density of states on the Ag (full line) and Se 
(dotted line) sites. The arrows indicate the energies of the 
states analysed in Fig. 3. 

In Fig. we show the electronic density of states 
(DOS) of a-AfeSe. The local densities of states (LDOS) 
shown in Fig. Ep allow us to identify the main features of 
the DOS. The Se(4s) states give rise to the lowest peak 
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at -12.5 eV, while the large peak below -4 eV originates 
from the Ag(4d) states. This prominent peak is super- 
imposed on a broader feature extending from -6.5 eV up 
to the Fermi level, which arises from the Se(4p) states. 
Even though the DOS is very small at the Fermi level, 
there is no real gap in the system. This is not surprising 
in view of the very small value of the experimental gap 
(less than 0.25 eV) and the well known tendency of DFT 
to give energy gaps that are too small. An examination 
of the band structure reveals that the states at the Fermi 
level are situated around the F point in the BZ. 

The LDOS reveals that the Se(4p) band is composed of 
two peaks situated at -5.5 eV and -1.5 eV. These peaks 
coincide with features in the LDOS on the silver atoms, 
which exhibits an asymmetric shoulder below its main 
peak and has a tail that extends from -3 eV up to the 
Fermi level. We believe this to be a signature of hy- 
bridization between the Ag(4(i) and the Se(4p) states. 
This hypothesis is clearly supported by examining den- 
sity plots arising from single bands. 




FIG. 3. Square of Kohn-Sham orbital from two bands with 
energies of -5.76 and -2.30 eV respectively. The orbitals are 
plotted on a plane passing through a Se (circle) and two neigh- 
boring Ag (square). 

In Fig. ^ we show the densities of two states, in a 
plane passing through a Se atom and two neighboring 
Ag atoms. The energies of the two states are indicated 
by arrows in Fig. |^. For both energies, the density is dis- 
tributed over both the Se and the Ag atoms. The wave- 
function clearly has a p character on the Se, whereas on 
the Ag sites the four visible lobes indicate a pronounced 
d-like character. 

An alternative way of looking at the bonding in a solid 
is to examine how the charge is redistributed with respect 
to neutral atoms. A careful analysis of this charge redis- 
tribution in Q;-Ag2Se shows that the region where the 
electron density has been enhanced is centered on the 
Se atoms and is broadly - almost spherically - spread 
between the Se and its nearest Ag atoms. This non- 
directional character of the bonds together with the high 
coordination of Se are compatible with ionic bonding. 

The picture that emerges is thus a valence band con- 
sisting of Se(4s), Se(4p) and Ag {4d) states with strong 



hybridization between the latter two. The empty states 
above the Fermi level can be seen as Ag{5s/p) states, pre- 
sumably hybridized quite strongly with Se(4p). With the 
Se(4p) band filled and the Ag(5s/p) states empty we thus 
arrive at a (partially) ionic model for the stoichiometric 
solid. 



IV. LIQUID ALLOYS 

A. Details of the simulations 

Our simulations of the Ag-Se liquid alloys have been 
performed on a system of 69 atoms in a cubic box with 
periodic boundary conditions. The wave-functions are 
expanded in plane waves with the same cutoff of 400 eV 
as before; this implies a basis set of ^ 28,000 plane waves 
for the whole system. We included only the F-point to 
sample the Brillouin zone, and we used a Fermi-smearing 
energy width of 0.2 eV, as for the solid. The Verlet algo- 
rithm was used to integrate the ionic equation of motion, 
with a time step of 3 fs. We have performed simulations 
at the temperature T ~ 1350 K for three concentrations 
of Agi-ajSe^:, namely: x = 0.33 (46 Ag atoms and 23 Se 
atoms), X = 0.42 (40 Ag atoms and 29 Se atoms) and 
X — 0.65 (24 Ag atoms and 45 Se atoms). The sim- 
ulations were performed at a density which is linearis 
interpolated between experimental values for £-Ag2SeE2l 
and i-Sdi. 

To initiate the simulations, we esploited the fact that 
an empirical pair-potential modelE3 has been developed 
for the stoichiometric Ag2Se system. This is a partially 
ionic model which reproduces the structure of the liquid 
reasonably when used in classical molecular dynamics - 
a direct comparison of the model with our AIMD simula- 
tions will be presented later. We began by making simu- 
lations with this empirical model, and we then switched 
over to AIMD and let the system equilibrate for a fur- 
ther 1 ps before collecting data over the next 2 ps. We 
reached the other Agi_a;Sea; compositions by replacing 
some of the Ag atoms by Se atoms and then equilibrat- 
ing for 1 ps at the new composition; production runs of 
1-2 ps were again performed in each case. 

B. Structural properties 

1. Structure factors 

The most direct and detailed way to compare our sim- 
ulations with experiment is through the static structure 
factors, which have been measured for the stoichiomet- 
ric composition Ag2Se. The neutron- weighted structure 
factor Sn{k) was measured some years agolij, and very re- 
cently the partial structure factors Sap{k) have also been 
measured, using the technique of isotope substitutionO. 
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The definition of Safsik) used here to ccmpare with ex- 
periment is that due to Faber and ZimarEj, according to 
which: 

Saf3ik) = {{p*^{k)pi3{k)) - 5ap)/^/c^j, (1) 

where Cq is the concentration of species a, Pa{k) rep- 
resents the Fourier transform of the number density of 
species a, and the angular brackets indicate the ther- 
mal average. With this definition, the neutron-weighted 
structure factor is given by: 

^ _ 'EallCgCisbabilSalsik) ^ ^ 

where ba is the coherent scattering length of species a. 
We usCriiie scattering lengths ^Ag — 5.98 fm and 6so = 
7.97 fJ^IM 

Comparisons of the structure factors given by simu- 
lation and experiment are shown in Fig. ^ Agreement 
between the two for S'„(fc) is very good, with all the main 
peaks having the correct position and height. The excel- 
lent agreement of the period, amplitude and phase of the 
oscillations at larger wavevectors is particularly gratify- 
ing. Our results also reproduce the small pre-peak at fc ~ 
1.7 A~^, though its magnitude is somewhat lower than 
the experimental value. 
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FIG. 4. (a) Comparison of the calculated total neutron 
weighted strueture factor Sn{k) of Ag2Se with the experi- 
mental resultaj. i-^ib) The measured partial structure factors 
Sa^(fc) of Ag2SeLJ compared with those extracted from our 
simulation. The theoretical curves are shown in thick lines, 
those from experiment with circles (the thin line is a smoothed 
version of the experimental curve). 

However, the more detailed comparisons allowed by the 
Safiik) structure factors reveal some discrepancies. The 
most noticeable is the absence of the pre-peak at fc ~ 
1.7 A-i in our SAg-Ag{k)- There are also substantial 
differences in the magnitudes of the peaks in S'Ag-So(fc) 



and Sse-Seik) at this wavevector. Since we found a dis- 
crepancy in Sn{k) at the same wavevector, the problem 
is unlikely to be due to the experimental analysis of the 
structure factor into partials. Our initial suspicion was 
that the small size of our system might be responsible for 
the discrepancy. However, we have tested this idea by 
doing simulations with the empirical model for different 
sizes of system, and these tests give no indication of sig- 
nificant differences between the 69-ion system and much 
larger systems. It is possible that fluctuations at the 
wavevector in question are very slow and that our simu- 
lations are not long enough to achieve adequate statistics 
for such long wavelengths. 

2. Pair correlation functions 

Our calculated partial radial distribution functions 
(RDF) 5Ag-Ag(?'), gAg-Se{r) and gsc-Scir) for the sto- 
ichiometric composition are reported in Fig. where we 
coroBare with the neutron diffraction results of Lague et 
al. lij, andnwith the predictions of the empirical interac- 
tion modelEHI mentioned earlier. 
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FIG. 5. Calculated partial radial distribution functions 
gapir) of AK2Se (full lines) compared with the results of 
experimentlij (d«tted lines) and the predictions of an empiri- 
cal ionic modelE3 (dash-dotted lines). 

The agreement of our simulated RDFs with experiment 
in the region of the first peak is excellent for Ag-Se and 
reasonably good for Ag-Ag, but in both cases there are 
quite noticeable differences at larger distances, as would 
be expected from the discrepancies of the corresponding 
Sap{k) at fc ~ 1.7 A^^. The agreement with experiment 
is somewhat less good for gsc-Sc{T), and there is no sign 
in the experimental results of the doubling of the first 
peak observed in our simulations, or of the weak dis- 
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tribution we find for r < 3 A. We return below to the 
meaning of the latter. If we are correct in attributing 
discrepancies in the structure factors at low k to inade- 
quate sampling, then discrepancies in the RDFs at large r 
are presumably due to the same effect. However, the dif- 
ferences of simulated and experimental 5se-Sc('') around 
the first peak must have a different explanation, and we 
do not at present know what this is. 

Our comparisons show that the empirical interaction 
model performs extremely well for gAgScif): but differs 
greatly from both the experimental and AIMD results 
for gAg-Ag('''). It seems certain that the model is se- 
riously mis-representing the Ag-Ag interaction at short 
distances. Interestingly, for gsc-Scir), the model agrees 
somewhat better with experiment than our AIMD results 
around the first peak. 
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FIG. 6. Partial radial distribution functions ga/iir) from 
simulated £-Agi_i,Sei, at concentrations a;=0.33 (full line), 
x—0.42 (dotted line) and x=0.65 (dot-dashed line). 

We show in Fig. || how the partial RDFs change with 
composition according to the AIMD simulations. The 
results show that increase of Se content causes dramatic 
changes in gsoSeir)- At stoichiometry, gsG-So(?') exhibits 
broad peaks at 3.99 and 4.83 A, and only a very weak 
tail below 3.0 A. For x = 0.42, instead of a tail, there is a 
short-distance peak at 2.35 A. The position of the main 
peak has shifted to 4.72 A and has decreased in magni- 
tude. The peak around 4.0 A seen at stoichiometry has 
merged with the main peak to give rise to a shoulder. At 
the last concentration, the short-distance peak is domi- 
nant and 5sc-So('') shows little structure beyond 4.0 A, 
with a low and broad second peak. 

The mechanism behind the changes observed in 
5Sc-Sg('') is the tendency of Se to bond to itself as x ex- 
ceeds 0.33, and it is no coincidence that the position of 
the short-distance peak is almost the same as the Se-Se 
covalent bond lengths_in crystallinq_and liquid Se, which 
are equal to 2.37 Ae3 and 2.34 All3 respectively. The 
nature of this peak can be probed by analysing the coor- 
dination of the selenium atoms. In general, we define the 
a-P coordination as the average number of neighbouring 



(3 atoms within a sphere of radius Tc around an a atom. 
The growth of the short-distance peak can be character- 
ized by the Se-Se coordination, which, for rc = 2.9 A, 
we calculate to be 0.1, 0.6 and 1.5 for x — 0.33, 0.42 
and 0.65 respectively. Thus the coordination grows as 
X increases, reflecting the tendency of Se atoms to form 
bonds. This will become even clearer when we examine 
the electronic structure of the liquid below. In compar- 
ison to the modifications observed in gse-Se{r), gAg-Sei^) 
is hardly affected by the change of composition. The 
magnitude of the main peak decreases slightly with in- 
creasing X and its position varies only a little: 2.65, 2.60 
and 2.64 A for a; = 0.33, 0.42 and 0.65 respectively. At 
stoichiometry, a very broad second peak may be distin- 
guished around '--^ 6 A, but it progressively disappears 
with increasing x, and at x = 0.65 gAgScif) is nearly 
constant for r > 4.0 A. A similar gradual loss of structure 
is observed for the Ag-Ag pair correlation function as the 
height of the first peak decreases with increasing Se con- 
centration. This is accompanied by a slight shift of the 
first peak of gAg-Agif) to smaller r. For x — 0.65 we note 
the presence of many oscillations in (7Ag-Ag('")- These may 
partly be due to the small number of silver atoms at this 
concentration and the resulting poorer statistics for the 
Ag-Ag pair correlation function. However, the absence 
of structure in gAg-Ag{f) at this concentration clearly in- 
dicates that the arrangement of the Ag atoms is very 
disordered. Furthermore, the low magnitude of the first 
peak suggests that the Ag atoms are in fast exchange with 
their first coordination shell, and that the Ag atoms are 
likely to be very mobile. The characteristic interatomic 
distances and coordination numbers deduced from the 
partial RDFs are summarized in Table together with 
those found in both the low temperature (a-Ag2Se) and 
the superionic (/3-Ag2Se) phases of Ag2Se. 

TABLE II. Inter-atomic distances and coordination num- 
bers in simulated l-Agi-xSe^ and in the crystalline phases 
Qf-Ag2Se and /3-Ag2Se. Va/s is the first neighbor distance be- 
tween species a and p. r'^p is the cut-off radius used for 
calculating the coordination number n„^. All distances are 
given in A. 



Phase 


X = 0.33 


X = 0.42 


X = 0.65 


a-Ag2Se 


/3-Ag2Se 


''Ag-Ag 


2.93 


2.84 


2.85 


2.93 


3.15 


'"Ag-Se 


2.65 


2.60 


2.64 


2.62 


2.70 


rSe-Sc 


3.99,4.83 


2.41,4.72 


2.35 


3.98,4.59 


4.30 


''Ag-Ag 


3.90 


3.93 


3.76 


3.68 


4.35 


'"Ag-So 


3.77 


3.65 


3.62 


2.90 


3.45 


'"Se-Se 


2.9,6.12 


2.90 


2.90 


5.14 


6.00 


lAg-Ag 


6.6 


3.4 


2.0 


6 


10 


flAg-Se 


3.8 


4.0 


5.1 


3.5 


4 


"■Sc-Sc 


0.1,13.6 


0.6 


1.5 


14 


14 
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It is interesting to notice that the distances and coordi- 
nation found for X = 0.33 are very cfose to those found in 
the fow temperature phase. In particular, the positions 
of the two superimposed peaks in 5sc-Sc('') are similar to 
the two nearest Se-Se distances found in the solid. This 
is surprising, since one would rather expect the melt to 
resemble the disordered superionic phase, in which the 
Ag ions diffuse through a b.c.c. Se lattice. 



3. Analysis of the Se structure 

The growth of the short-distance peak in gsc-Scif) 
with increasing Se content is due to the formation of Se 
clusters for x > 0.33. The simplest way to see this is 
by studying 'ball-and-stick' pictures for typical configu- 
rations. Fig. 1^ displays examples of such pictures for the 
three compositions. In constructing these pictures, pairs 
of Se atoms are joined by sticks if their separation is less 
than 2.9 A, this distance being chosen because it is close 
to the first minimum of gsc-Scir) for x = 0.42 and 0.65. 
At stoichiometry, very few Se clusters are present, with 
only the occasional dimer (Fig. ^), and the even more 



(a) 



(b) 




FIG. 7. Snapshots of typical configurations of l-Pi-g\-x^&x 
at concentrations (a) 1=0.33, (b) x=QA2 and (c) x=0.65. Sil- 
ver atoms are shown as black spheres, selenium atoms as gray 
spheres. Bonds are drawn between Se atoms with separation 
< 2.9 A. Bonds to atoms in neighboring cells are represented 
with two-colored sticks. 

occasional trimer. It is the presence of these small clus- 
ters that gives rise to the weak tail in S'Se-Se('*) at low 
r. For X = 0.42, the Se atoms bond not only into dimers 
and trimers but also into larger Se„ clusters (see Fig. 
^p). At this composition, some ~ 48 % of the Se atoms 
are bonded, with 76 % of these one-fold and 22 % two- 
fold coordinated, with the remaining 2 % having higher 
coordination. 

TABLE III. Distribution of coordination numbers nsc-Sc 
at a;=0.33, 0.42 and 0.65. We show the average percentage of 
atoms with coordination Nc. 



Ae 


X = 0.33 


X = 0.42 


X = 0.65 





93 


52 


9 


1 


6 


36 


46 


2 


<1 


10 


38 


3 




<2 


6 


>4 






<1 



The dominance of one-fold and two-fold coordination 
means that most of the clusters are either dimers or 
Se„ chains. At the composition x = 0.65, most of the 
Se atoms are in clusters, with an average of 7 % be- 
ing isolated (see Fig. ^). The proportions of bonded Se 
with one-fold and two-fold coordination are now 35 % 
and 40 %, with a rather significant percentage (15 %) of 
three-fold coordinated atoms. In Table III, we summa- 



rize the distribution of first-neighbor Se-Se coordination 
numbers, for all three concentrations. 

To understand the structure of the Se clusters in more 
detail, we have studied their size distribution. In Fig. H, 
we display the average numbers of clusters present in our 
sample versus the number of Se atoms n in the cluster. 
For all three concentrations, dimers are the most com- 
mon clusters. The number of clusters drops rapidly with 
increasing n. For x ~ 0.42, the prevalent clusters contain 
up to 5 atoms. Only a very small number of clusters with 
n > 5 was found, the largest being n = 9. When the Se 
content increases to x = 0.65, the cluster-size distribu- 
tion broadens substantially. The largest cluster present 
in our sample was n = 28, but all clusters with n < 10 
are very frequent. 




FIG. 

^-Agi-. 
line) . 



:Se= 



8 12 

Cluster size 

Average number of Se clusters of size n in 
at concentrations 2;=0.33, 0.42 and 0.65 (full 



An analysis of the structure of the clusters reveals 
that there is a simple relationship between the number 
of bonds within a cluster Ni, and the number of atoms 
composing that cluster n. We found that for over 95 % 
of the clusters we have A^b = ^ 1, at all three con- 
centrations. The significance of this is that the relation 
Nh = n — 1 indicates a tree-like topology without closed 
loops. Almost all the exceptions to this are clusters with 
Ni, = n, which indicates the existence of a single closed 
loop. We note that, for large clusters, loops defined in 
this way may occur simply as an artifact of the periodic 
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boundary conditions. In any case, the conclusion is clear: 
loops are rare. 



C. Dynamical properties 

We have studied the diffusion of atoms in the liquid 
by calculating the time-dependent mean square displace- 
ment (MSD), defined in the usual way for species a as: 




(3) 



where the sum goes over all atoms of species a, to is 
an arbitrary time origin, and the angular brackets denote 
a thermal average, or equivalently an average over time 
origins. For diffusing systems, the MSD is linear in t for 
large \t\, and the slope is proportional to the self-diffusion 
coefficient Da of species a: 



(Araitf) ^ 6Da\t\ + 
where is a constant. 



(4) 
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FIG. 9. Mean square displacement (Ara(t)^) of Ag (top 
panel) and Se (lower panel) in ^-Agi-^Sea; for a;=0.33 (full 
lines), 0.42 (dotted lines) and 0.65 (dash-dotted lines). 

We have calculated the MSD for Ag and Se for the 
three compositions, by averaging over the atoms of each 
species and over time origins. Our calculated MSDs are 
graphed in Fig. 0, and the estimated diffusion coeffi- 



cients are reported in Table IV. In all normal liquids, 
the asymptotic linear behavior of (Ar„(t)^) is attained 
after only a few tenths of a ps, and this is the behav- 
ior we find. The lack of straightness of the curves at 
long times is simply an effect of statistical averaging (the 
ideal asymptotic form shown in Eqn. assumes perfect 
averaging over a simulation of infinite duration). The 



slopes used to obtain the values of Da were estimated by 
a least-square fitting of a line to the data. 

Our results show three things: the diffusion coeffi- 
cients are similar to those of other normal liquids, be- 
ing in the range 10~^ - 10~^ cm^s~^; Ag diffuses faster 
than Se; and the diffusion of both species increases with 
Se content. Unfortunately, there appear to be no ex- 
perimental data for the diffusion coefficients. However, 
Z?Ag lias been measured in the superionic phase of the 
/3-Ag2Se crystalHEj. This is relevant, because it is com- 
monly found that the diffusion coefficients of the mobile 
species in superionic conductors change only a little on 
melting. The measured value of i^Ag in superionic /3- 
Ag2Se at the melting point is 6 x 10~^ cm^s~^E3, which 
is rather close to our calculated value of 6.4x 10"^ cm^s~^ 
for the stoichiometric liquid. It is also relevant to note 
that the empirical interaction modelEj reproduces the ex- 
perimental values of -DAg in /3-Ag2Se quite well. This 
model predicts the values I?Ag ~ 10.6 x 10~^ cm^s 
Dsc = 3.4 X lO-^cm^s-i for ^AgzSe at - 1380 KB 
which are quite close to our AIMD values. However, 
we believe that too much weight should not be put on 
this comparison, because of the rather poor results for 
gAg-Agir) given by the empirical model. 

Even though experimental values for the Da are not 
available, there have been measurements of the ionic con- 
ductivity (7i, i.e. the electrical conductivity measured un- 
der conditions such that the flow of electrons is blocked. 
A rough check against our calculated Da can be made 
using the approximate Nernst-Einstein relation: 



CFi = {kBT) ^y^^PajzaefDq 



(5) 



where pa is the bulk number density of species a and Za 
is its ionic charge. This relation is based on the assump- 
tion that ions diffuse independently of one another. The 
check can only be rough since it is not clear exactly what 
values should be assumed for Za- However, if we assume 
full ionicity (zAg — -1-1, ^Sc = —2), then our calculated 
values of Da for the stoichiometric composition lead to 
a predicted ai of 4.9 ri~^cm~^, whiC|h_is-|rather close to 
the measured value of 5.2 f7^^cm^^Eli£3'Ea. 

TABLE IV. Self-diffusion coefficients Da of Ag and Se in 
^-Agi-^Se^, for a;=0.33, 0.42 and 0.65 in units of 10"^ cm^s"^ 



X 


0.33 


0.42 


0.65 




6.4 


11.2 


13.1 


Dsc 


2.7 


4.2 


6.3 
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Our finding that DAg and Dse increase with Se con- 
tent is surprising at first sight, since we have seen that 
Se forms large clusters when x > 0.3. As noted in the In- 
troduction, the chain-like nature of pure £-Se is also well 
established. It would be expected that the formation of 
covalently bonded clusters would make it more difficult 
for Se atoms to diffuse, and that the chain-like network 
would also hinder the diffusion of Ag. However, this de- 
pends very much on the stability of the chains, in other 
words on the lifetime of the covalent bonds. 




FIG. 



10. 



1 

Life time (ps) 

Distribution of Se-Se 
for x=0.33, 0.42 and 0.65. 



bond lifetimes in 



In order to pursue this question further, we have ana- 
lyzed our simulations to estimate the Se-Se bond lifetime. 
To define the lifetime of a bond between two atoms, we 
regard the bond as being made when the distance be- 
tween the atoms becomes less than a creation distance 
ri ; it is broken when the distance becomes greater than a 
somewhat larger cut-off r2 ■ We decided to adopt this def- 
inition involving two cutoff radii because there are large 
fluctuations in the distance between bonded atoms, and 
we do not wish to count a bond as broken simply because 
the distance momentarily exceeds the creation distance 
ri. In fact, our conclusions would not be significantly 
affected if we put ri = r2 . The distribution of bond life- 
times obtained for n = 2.9 A and r2 — 3.2 A is shown for 
all three compositions in Fig. 



10, The distributions are 



very broad, with lifetimes ranging from less than 0.1 ps 
up to 2.6 ps. For the majority of bonds, the lifetime is less 
than 1.0 ps. For x = 0.42 and 0.65, the average lifetimes 
are 0.52 and 0.50 ps respectively. (The number of Se-Se 
bonds in the stoichiometric case was too small to give 
a statistically meaningful lifetime.) These rather short 
bond lifetimes make it less surprising that the diffusion 
coefficients increase with increasing Se content. 



D. Electronic properties and bonding 

In this section we will examine how the structural 
changes are linked to the electronic properties of the liq- 
uid alloy. The calculated electronic DOS for the three 
compositions are shown in Fig. |ll|a. The states can be 
classified as: Se(4s) states around -12 eV, Ag(4(i) states 
at -4 eV and Se(4p) states in the region above -7 eV. 
The identification of these features is made straightfor- 
ward by the examination of the LDOS on the Ag and Se 
atoms, which are displayed in Fig. [TTIb. 

The DOS and LDOS for £-Ag2Se are very similar to 
those found for Q;-Ag2Se (see Sec. III). In particular, the 
hybridization of Se(4p) states with Ag{Ad) states is also 
found in the liquid, as is clear from the asymmetry of the 
Ag(4d) band and the dumbbell-shaped Se(4p) band. We 
note however an increase in the DOS in the region of the 
Fermi energy on going from solid to liquid. 
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FIG. 11. Density of states (left panel) and local densities 
of states (right panel) for Ag (chain curve) and Se (dotted 
curve) from simulations of Agi-ajSe^; at j: = 0.33, 0.42 and 
0.65. For clarity, the scale used for the Se LDOS is four times 
that used for the Ag LDOS. The vertical dotted line marks 
the Fermi energy. 

Even though the DOS for the three concentrations are 
quite similar, we observe significant modifications as the 
Se concentration is changed. First, the height of the 
sharp Ag(4d) peak decreases with respect to the higher 
peak of the Se(4p) band which remains the same. Sec- 
ond, the Se(4s) band which is narrow and has a single 
peak at stoichiometry becomes substantially broader and 
splits into two peaks as x is increased. Finally, the DOS 
at the Fermi level {E-p ) increases with x and the pseudo- 
gap present at a:=0.33 becomes shallower. 

The LDOS reveal the nature of these changes. We see 
that the shape of the Ag LDOS is essentially identical for 
all three concentrations. In particular, the asymmetry of 
the main peak is present for all x, which means that the 
Ag(4(i)-Se(4p) hybridization is unaffected by changes in 
the concentration. The rigidity of the Ag LDOS thus im- 
plies that the modifications in the DOS are mainly due 
to the changes in the LDOS of Se. There are indeed 
large changes in the Se LDOS as x changes from 0.33 to 
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0.65: the two peaks of the p band are smeared out and a 
third peak appears above the Fermi level. The first and 
third peaks are situated symmetrically about the mid- 
dle peak and are approximately of the same magnitude. 
The structure of the Se LDOS at x = 0.42 may simply 
be described as intermediate between these two cases, as 
can be seen from the Se(4s) states. For all three con- 
centrations, we see that the states at E-p arise from the 
Se(4p) band. Thus the increase of the DOS at the Fermi 
level is clearly due to the broadening of the second Se(4p) 
peak, and to the presence of the third peak which both 
effectively reduce the depth of the pseudo-gap. 
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FIG. 12. Decomposition of the Se LDOS into contributions 
from isolated (OF), one-fold (IF), two-fold (2F) and three-fold 
(3F) coordinated Se atoms. 

All these changes are closely related to the modifica- 
tions of the structure. Our interpretation of the changes 
is as follows. The presence of a pseudo-gap and the low 
value of the Se LDOS above the Fermi level at the stoi- 
chiometric composition suggest that all the Se(4p) states 
are filled, so that there are 6 p-electrons per Se atom. 
This corresponds to an ionic model in which it is ener- 
getically favorable for the Ag(5s) electrons to transfer 
to empty Se(4p) states, since the atomic Ag(5s) state is 
higher in energy than the Se(4p) states. This corresponds 
to a charge transfer from the Ag to the Se ions, thus lead- 
ing to Ag+ ions and Se^^ ions. This naturally leads one 
to interpret the states above E-p as being Ag(5s-p) states, 
so that the picture here is essentially the same as for the 
solid, where we have already seen that an ionic model is 
appropriate. As x increases beyond 0.33, the number of 
occupied Ap states per Se atom falls below 3 due to the 
creation of Se-Se bonds. The formation of these bonds 
is clearly illustrated by the appearance of the third peak 
above the Fermi level, which corresponds to the anti- 
bonding states, the first and second peak corresponding 
to bonding and lone-pair states, respectively. The anti- 



bonding states are hybridized with Ag(5s-p) states, and 
only become clearly visible at a; = 0.65. The formation 
of Se-Se bonds is also reflected by the splitting of the 
Se(4s) which arises from the formation of bonding and 
anti-bonding Se(4s) combinations in the clusters. 

To show more clearly the effect of the formation of 
Se-Se bonds on the Se LDOS, we have decomposed this 
LDOS for X = 0.65 into contributions from isolated, one- 
fold, two-fold and three-fold coordinated Se atoms, as 
shown in Fig. |2| (we denote these different coordinations 
as OF, IF, 2F and 3F). We note the strong resemblance 
of the LDOS for the isolated atoms with the Se LDOS for 
the stoichiometric composition, with the Fermi close to 
minimum of the LDOS. For 2F atoms, the Se(4s) band 
has split symmetrically into bonding and anti-bonding 
states, while the Se(4p) now consists of three peaks, 
which we have identified above as bonding, non-bonding 
(lone pair) and anti-bonding. The Fermi energy in this 
case falls roughly between non-bonding and anti-bondiac' 
states, as we would expect in a chain-like structureE3. 
The LDOS for IF atoms is intermediate between the OF 
and 2F forms; the Se(4p) bonding and anti-bonding peaks 
are weaker than for the 2F atoms, since only one bond 
is involved. Note that if Ep fell exactly at the top of the 
non-bonding peak the number of electrons in the occu- 
pied Se(4p) states would give the OF, IF and 2F atoms 
charges of -2, -1 and respectively. 




FIG. 13. Square of Kohn-Sham orbital from two individual 
bands with energies of -5.3 (a) and -1.0 eV (b) respectively. 
The orbitals are plotted on a plane passing through a Se (cir- 
cle) and two neighboring Ag (square). The hybridized nature 
of the states is clearly visible, with a p and d character of the 
wave functions on the Se and Ag atoms respectively. The hy- 
bridization has a bonding character for the lowest states and 
is anti-bonding for the highest states. 

Several features of the DOS and LDOS can be con- 
firmed by a study of the electron density. A useful way 
to do this in the liquid is to examine the density in a plane 
passing through three neighbouring atoms. The presence 
of hybridization of the Se(4p) states with Ag(4(i) states 
is illustrated in Fig. ^ which shows the density of one 
state from the bottom and another from the top of the 
Se(4p)/Ag(4(i) band at -5.3 eV and -1.0 eV, respectively. 
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The very close resemblance of these density distributions 
to those found in the crystal (see Fig. ^ shows that this 
aspect of the electronic structure persists essentially un- 
changed in the liquid state. In the ionic model we pro- 
pose the states in the Ag LDOS above the the Fermi level 
arise from the empty silver s states. This is confirmed 
by the charge density of a state at 2 eV above the Fermi 
level shown in Fig. |lj. The wave function has a clear s 
character on the Ag atom and is slightly hybridized with 
Se(4p) states. 




X 

FIG. 14. Square of Kohn-Sham orbital of a state from the 
conduction band with an energy of 3.5 eV. The orbitals is 
plotted on a plane passing through a Ag (square) and two 
neighboring Se (circle). This state is mainly localised on the 
Ag atoms and has a clear s character. This supports the ionic 
model in which the silver atoms have lost their s electrons to 
become Ag"*" ions. 

We have calculated the electronic conductivity jO. 
in our simulations by the Kubo-Greenwood formulaEll 
as has b(pti— . dope in several previous ab initio 
calculationsHcjSo. Our results are shown in Table 
where they are compared with recent experimentaLmea- 
surements made over the whole composition rangeLrLl. We 
have included in the table a calculated result for the con- 
ductivity of -^-Se taken from ourxiwn unpublished AIMD 
simulation performed at 1370 KL2I, since this will be rel- 
evant to the discussion. 



TABLE V. Calculated conductivity (units of Sl^^cm^^) 
using the Kubo-Greenwood formula compared with experi- 
mental dataQ at different compositions x. 



X 


Theory 


Expt. 


0.33 


480 


410 


0.42 


1500 


300 


0.65 


1250 


490 


1.00 


400^^ 


lb 



The conductivity a = 480ri~"'^cm~^ found at stoi- 
chiometry is in satisfactory agreement with the exper- 
imental data. However, as x increases, the calculated 
a first increases dramatically to ~1500 il^^cm"^ at 
X — 0.42 and then decreases slightly to ~ 1250 fi^^cm^^ 
at a; = 0.65. These values are 3-4 times higher than the 
experimental values. For pure ^-Se the calculated a is 
too large by at least two orders of magnitude. We believe 
that these large discrepancies with experiment arise from 
the approximations implicit in the Kubo-Greenwood ap- 
proach. There are two major deficiencies in the Kubo- 
Greenwood formula as applied here. First, it completely 
ignores electron-electron collisions, and includes only the 
effect of scattering of the electrons by the ions. Second, 
and much more serious, it works with the Kohn-Sham 
energies, rather than with quasi-particle energies. For 
metallic systems having free-electron-like DOS— the er- 
rors_are probably small, and in cases like f-SiE3 and l- 
GaE3 the method has been shown to work well. But in 
systems having a pseudogap, like the liquid alloys consid- 
ered here, the use of the (artificial) Kohn-Sham energies 
is very likely to cause serious errors. The case of pure l- 
Se is particularly instructpe, since experimentally there 
is a band-gap of ^^0.7 eVtJ under the conditions of in- 
terest, whereas our AIMD simulations show no gap, but 
only a deep minimum in the DOS. It is important to be 
clear that this does not point to any error in the DFT 
calculations as far as liquid structure and dynamics are 
concerned. What is at issue here is the (unsurprising) 
consequence of trying to use Kohn-Sham energies in an 
unjustified way. Similar large overestimates of the elec- 
tronic conductivity havCr-been reported in recent AIMD 
simulations on ^-Mg3Bi2Ej. However, it is not yet clear 
under what circumstances one should expect these large 
errors, and it is worth noting that the conductivity calcu- 
lated in our very recent AIMD simulations on the £-Ga- 
Se systemcj agreed quite well with experiment in spite of 
their low values and the existence of a deep pseudogap. 

V. ATOMIC ORDERING AND ELECTRONIC 
PROPERTIES 

We have remarked in a number of places on the ionic 
nature of the solid and liquid Ag-Se system, and the im- 
plications of this for understanding the structure. We 
now try to bring these remarks together into a coher- 
ent picture. The ionic nature of solid Q!-Ag2Se is clear 
from an analysis of the DOS, and we have seen that the 
electronic structure changes little an melting, so that l- 
Ag2Se can be regarded as consisting approximatively of 
Ag+ and Se^~ ions. In the more Se-rich liquid we have 
used the DOS and LDOS to argue that isolated Se ions 
are in the Se^~ state, while the chains consist of neutral 
Se atoms terminated by Se~ dangling bonds. This means 
that the hyperstoichiometric liquid consists of Ag+, Se'^" 
and (Se„)^~ complexes. This picture can be tested by us- 
ing simple charge-balance arguments to predict the aver- 
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age number of Se-Se bonds in the liquid at any compo- 
sition, and hence the coordination number of the short 
distance peak in (7Sc-Sc('')- If the system consisted en- 
tirely of Ag+ and Sc^^ ions, the net charge on a system 
of iVAg and A^'gc ions would be iVAg — 2A^go. If the clusters 
are all (Se„)^" then the formation of every Se-Se bond 
reduces the charge by two units. For electroneutrality, 
the number of bonds must therefore be iVgo — ^-^Ag, and 
the coordination number of the short distance peak in 
5Sc-Sc('') must be 2 — Nj^g/Ngc — 3—l/x. This formula 
gives predicted values of 0, 0.62 and 1.46 for the Se-Se 
coordination number at the compositions x = 0.33, 0.42 
and 0.65, which are quiet cl ose to t he values 0.1, 0.6 and 
1.5 that we reported in Sec. IVB 2. This agreement con- 
firms the picture we propose. 

We can interpret the composition dependence of the 
experimental electrical conductivity in the light of this 
picture. For pure i-Ag, the DOS is expected to be rather 
free-electron-like near Ep. As Se is added, electrons are 
drained from the Ag(5s/p) band into the Se{4p) states 
until, at the stoichiometric composition, Ep is in a pseu- 
dogap at the top of the Se(4p) band. If the DOS remained 
rigid beyond this point, Ep would move into a region 
of higher DOS, and a would increase strongly. Instead 
of this, Se-Se bonding starts to occur, and this creates 
a pseudogap between Se non-bonding and antibonding 
states. It is this that maintains a at low values. Finally, 
as pure £-Se is approached, this pseudogap widens, and a 
falls to even lower values. As we have seen in the previ- 
ous section, our numerical values of a do not reproduce 
this behavior quantitatively in the hyperstoichiometric 
region. 

If this picture is correct, it should also be rather gen- 
eral. It is worth noting that chalcogen pairing has been 
observed by diffraction exaeriment ip-,some equiatomic 
liquid mixtures, e.g. CuSeO and KTeC3. The conductiv- 
ity of £-Cu-Se has been measured over a range of compo- 
sition and behaves in a similar way to that found in our 
simulations, showing a strong minimum at the composi- 
tion Cu2Se and then passing through a maximum in the 
region of CuSe before falling to the low values associated 
with £-Se. It would clearly be interesting to investigate 
these and other related systems by AIMD simulation. 



tent exceeds x = 0.33, and this consists of the covalent 
bonding of Se atoms to form clusters. These clusters 
are mainly chain-like, but ed x — 0.65 there is a sig- 
nificant fraction of three-fold coordinated Se atoms; the 
concentration of rings is extremely small. The equilib- 
rium fraction of Se present in the form of clusters can be 
understood on a simple charge-balance argument based 
on an ionic interpretation. In spite of the Se clustering, 
the diffusion coefficients of both Ag and Se increase with 
Se content. This is not inconsistent with the existence of 
the clusters, because the Se-Se bonds turn out to have 
very short lifetimes of less than 1 ps. The electronic den- 
sity of states of liquid Ag2Se closely resembles that of the 
solid. With increasing Se content, the main changes in 
the density of states are the splitting of the Se(4s) into 
bonding and anti-bonding states and the formation of 
Se(4p) bonding and anti-bonding states associated with 
Se chains. The calculated electronic conductivity is in 
fair agreement with experimental values for the stoichio- 
metric liquid, but is too large by a factor of 3 - 4 for 
higher Se contents. We have suggested that this error 
comes partly from the use of Kohn-Sham energies rather 
than quasiparticle energies in the conductivity calcula- 
tions. 
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VI. CONCLUSIONS 

The AIMD simulations we have presented lead us to 
the following conclusions. The reliability and realism 
of the DFT-pseudopotential methods we have used are 
demonstrated by the satisfactory predictions of the equi- 
librium structure of the a-Ag^Se crystal a the quite close 
agreement with neutron diffraction results for the par- 
tial structure factors and RDFs. Nevertheless, there are 
significant discrepancies, which need further investiga- 
tion. The simulations have shown that a large structural 
change in the Agi_a;Sea; liquid begins when the Se con- 
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